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Abstract 

The present Standard Model prediction for muon g — 2 is reviewed. Em- 
phasis is put in discussing the main hadronic uncertainties. 
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1 Introduction 



Recently, the Muon g — 2 Collaboration from the E821 experiment at Brookhaven 
National Lab (BNL) reported a new result for the muon g — 2 with an un- 
certainty more than five times smaller than the last CERN experiment [0]. The 
E821 result when combined with the previous experiments produce the present 
world average 

= ^il-^ = (11 659 202.3 ± 15.1) ■ 10"^° . (1) 

The expected impressive final goal of E821 is to achieve an experimental uncer- 
tainty in of the order of 4 ■ 10"^''. 

Accompanying this great experimental performance a lot of effort has been 
put in the theoretical side to get a Standard Model prediction for this quantity 
since the pioneering work of Schwinger In the next Sections I will review 
the present Standard Model prediction putting emphasis in discussing the main 
hadronic uncertainties which at present are the dominant. 

2 Quantum Electro- Weak-Dynamics Contribu- 
tions 

2.1 QED Contribution 

The QED contribution to muon g — 2 is known to 0{a^) 

\ (^^) + 0.765 857 376(27) + 24.050 508 98(44) 

126.07(41) + 930(170) (-)^ + " " " (2) 

Higher orders are negligible compared with the experimental uncertainty. The 
original references leading to this result can be found in ^, ^. The first, 
second, and third orders are known fully analytically. The fourth order contains 
contributions which are only known numerically 0. The fifth order is a numerical 
estimate of the dominant diagrams enough for the BNL expected uncertainty. 
Using the present world average value = 137.035 999 76(50) 0, one gets 

^QED ^ (^^^ ggg g ^ Q_3^| . ^Q-IO _ (-3^ 

which gives the bulk of the experimental value of a^. 

^Recently, T. Kinoshita has found a numerical error in the fourth order contribution which 
final result is not public yet. However, the total Standard Model uncertainty is not upset by it. 
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2.2 Electroweak Contribution 



This contribution can be written as 



5 Gprnl 
38^2 



l + -(l-4siny%(-97.0±8.8) (^) + 



(4) 



where the first term is the one-loop contribution and the second term is the 
O^Gpm'j^a) two- loop contribution Using [Q, one gets 



;i9.4808--- + (-4.4±0.4)) ■ 10 



-10 



;i5.1±0.4) • 10 



-10 



(5) 



for the one-loop plus two-loop EW contributions. The uncertainty in the two-loop 
result takes into account uncertainties in the Higgs mass, quark two-loop effects 
and and higher corrections. The leading ln(M^/m^) logs can be re-summed 
to all orders in a in [|1^ and the result enlarge slightly the two-loop number but 
it is well within the quoted uncertainty, so that we keep the full two- loop result. 



3 Hadronic Contributions 



3.1 Hadronic Light-by-Light Contribution 

This is the contribution of a hadronic Green's function with four legs coupled 
to electromagnetic two-quark currents. This four-point function is attached in 
all possible ways with three of its legs to the muon line. This contribution is of 
0{a^) and cannot be related to any measured quantity and we have to rely in 
our ability of treating the strong interactions at all energies. 

There are always two topologies to a full four-point function; namely, first, 
two three-point form factors joined with full propagators, these we call three- 
point-like contributions, and second, the pure four-point form factor which we 
refer to as the four-point-like contribution. In all cases the three vector legs are 
joined to the muon line through full vector two-point functions. We calculate 
the leading 0{Nc) contributions to the full four-point function as well as NLO 
in 1/Nc corrections which are saturated by four-point-like charged pion and kaon 
loops. Being of different order in 1/iVc, they do not have to match the leading 
order in 1/Nc contributions as the quark- loop contribution [0, |12 



Based in the 1/Nc analysis of |]TT|, there have been two full calculations of 

T3[| . Here, I pay more attention to 



this contribution in n^A and 



12 and a full 



comparison with |TE[ and references to previous work can be found in |T2 



14 



Recently, there have also calculations of the pseudo-scalar exchange in 
and of the pseudo-scalar and scalar exchanges in DTSf. In [jl6l, the coefficient of 
the leading divergent logarithm has been obtained analytically for the point-like 



Wess-Zumino-Witten vr 77 form factor. This result has been confirmed in |17 . 
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The results we get are |T2| p| 



a^-^-^ = ((2.1 ± 0.3) + (-0.68 ±0.2) + (8.5 ±1.3) + (0.25 ±0.1)) • 10"^° 
± (-1.9 ± 0.5) ■ 10"^° = ((3.8 ± 0.3) ± (4.4 ± 2.1)) ■ 10"^° 
= (8.3 ± 3.2) • 10"^°, (6) 

where the first term is the four-point- hke contribution, the second, third, and 
fourth terms are the three types of three-point-hke contributions, namely, when 
the propagators are scalar, pseudo-scalar, or axial-vector. The first term in the 
second line is the NLO in 1/Nc contribution. The second split gives the contri- 
butions below 0.5 GeV and above 0.5 GeV, respectively. 

is 



With the same split as (pD, the result obtained in [13 



a^-^~^ = ((0.97 ± 1.1) ± (0.0 ±0.0) ± (8.27 ±0.60) ± (0.17 ±0.1)) ■ 10"^° 
± (-0.45 ± 0.8)) ■ 10"^° = (8.96 ± 1.54) • 10^^°. (7) 

Two main features appear from both results. The dominant contribution by far 
is the pseudo-scalar exchange and second, there is a very good agreement in the 
contribution from the pseudo-scalar exchange while for the rest of the 0{Nc) 
contributions the disagreement is pretty small. This disagreement is larger in the 
NLO in 1/Nc contribution. In fact both contributions cancel out very much and 
the full result is in quite good agreement. Further scrutiny of this cancellation in 
a as much as possible model independent way is needed. 

For the low energy part (below 0.5 GeV) of the four-point- like and scalar 
exchange three-point-like contributions -they are related by Ward identities- we 
use the ENJL model [|l^]. For the higher energy part of the four-point-like con- 
tribution we use a heavy quark which mass acts as an IR cut-off . The ENJL 
model has not free parameters, they are fixed to low-energy vrvr data. In the ENJL 
model, full vector two-point functions have an nonphysical behavior O, at 
intermediate energies which we corrected in this work. 

The pseudo-scalar three-point-like contribution is dominated by the 7r° ex- 



change with non-negligible contributions of the rj and rj' exchanges [12|. We used 
a variety of 71^7*7* form factors fulfilling all the known constraints. Fortunately, 
this form factor is very constrained by the Ua(1) anomaly which gives its normal- 
ization at the origin and by p° — *• 7r°7 which gives the slope at the origin. There 
are also data on 7r° 7*7 between 0.5 GeV and 3.3 GeV [^, All these 
constraints make the model dependence small and it is also the reason of the 
good agreement of (^ and (|^). Nevertheless, the uncertainty in this contribution 



can be reduced using data on e~^e — > 7r°e"'"e at intermediate energies |2^. Data 
below 0.5 GeV on vr'' 7*7 can also help. 

•^Notice that we correct for an errata in the sign of the pseudo-scalar and axial-vector ex- 
change contributions. 
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The authors of were able to obtain analytical formulas for the pseudo- 
scalar exchange for a general class of 71^7*7* form factors that fulfill the OPE 
and large Nc QCD constraints, and which are compatible with the data. Their 
result for this contribution (8.3 ± 1.2) ■ 10"^'' agrees very well with both and 



1^ after correcting for the sign mistake. 

For the NLO in l/N^ contribution we cannot use the ENJL model. We 
saturate them with charged pion and kaon loops coupled to photons and need 
P"'"P~7*7* and P+P~7* vertices for which we take complete VMD which works 
very well for one-photon couplings at all energies. For two-photons there is no 
data beyond tt+tt" 77. In [|T^, a HGS symmetry model was used and the dis- 



crepancy with [jT2| already at low-energy is large but being of CHPT order there 
is no data to disentangle it. It is possible that information on e+e^ 7r+7r^e+e^ 
could help to eliminate this large model dependence. 

The uncertainty in (j^) is obtained by adding linearly the individual errors 
plus 0.8 • 10~^° added linearly to take into account the discrepancy between our 
result and the one in (|^). The uncertainty in (|^ is obtained adding quadratically 
their individual errors. 

I take the average of @ and 

^L-b-L ^ ^g_g ^ 3 . IQ-W (g) 

as the present value for this hadronic light-by- light contribution. The uncertainty 
is the one from (j^) since already takes into account the discrepancy between both 
results as explained above and is more realistic. Improving both in the '7r''7*7* 
and 7r"'"7r~7*7* vertices could decrease this uncertainty to around 2 ■ 10~^°. 



3.2 Hadronic Vacuum Polarization Contribution 

This contribution starts at order and it is the one with the largest uncer- 
tainty at present. It is a hadronic Green's function with two legs coupled to 
electromagnetic two-quark currents and attached to the muon line. 
The 0{a^) contribution can be written as |]23[| , 



,(2)hvp 



4m1 



dt Kit) 



cr 



(0) 



hadrons) (t) 



a 



(0)| 



(9) 



where K{t) is a know function of t. There have been many calculations with 
increasing accuracy due to better data and theoretical input, see EM for recent 
calculations and 



25 



for a critical review. References to previous work can be 
found in ^, ^ Since its uncertainty can be reduced systematically with 
accurate data one should consider this contribution mainly of experimental origin. 

The most recent e~^e~ data |^ from BES-II at Beijing, and SND and CMD-2 
at Novosibirsk has been used in [B^ 



hvp 



(697.4 ±10.5) ■ 10 



'10 



(10) 
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Also just using e~^e data gets 

^e+e-hvp _ (698.2 ±9.7) ■ 10"^°. 
These results give the present average for this contribution 



a 



e+e hvp 



(697.8 ± 10.5) • 10 



^10 



(11) 



(12) 



Further reduction of the uncertainty in the e~^e data below the tau mass to the 
order or below 1% is expected from VEPP-2M (CMD-2, SND) at Novosibirsk, 
DA$NE (KLOE) at Frascati, BEPC (BES) at Beijing and other low energy 
facilities. In the theory side, the possibility of measuring the pion form factor at 



DA$NE and discussion of the relevant QED corrections was presented in [ p9|| . 
The complete 0{a) QED initial state, final state, and initial- final state radiation 
corrections to e^e~ tt+tt" has been recently presented |30|. 

Using CVC, one can relate the r~ — 7r~7r°i/^ with 1=1 vector channel 
data to the e^e~ tt~^7t~ data. The very precise tau hadronic decay data |^ 
from ALEPH and OPAL at CERN and CLEO at CESR when supplemented 
with the SU(2) breaking corrections can help in increasing the accuracy of the 
e'^e" — TT~^'K~ data which gives 70% of the contribution to and 80% of its 
uncertainty when integrated in (P) from 4m^ to t = 0.8 GeV. This program was 
started in and continued in where pQCD was also pushed down up to 1.8 
GeV helping to reduce the final uncertainty sizeably. This use has been recently 



confirmed by the BES-II data which give compatible results p3[] . 

Adding the rest of the hadronic vacuum polarization contributions and the 
known isospin breaking corrections, gets 



r hvp 



(692.4 ±6.2) ■ 10 



-10 



(13) 



The leading isospin breaking corrections studied in |^ agree quite well within 
errors with the same corrections applied in |3^. Numerically, the isospin correc- 
tions applied in p8| to the tau data are very similar to the ones in p3|. The 



result from tau data in PHI is 



a^^^^P = (695.2 ± 6.4) ■ 10" 



10 



(14) 



From these results the present average for the from tau data reads 

al^^P = (693.8 ± 6.4) ■ 10"^° . (15) 



The recent calculation is also compatible within errors though no isospin 
breaking effects are included. For the combined final result of (|T^ and ( jTSP I 
take the weighted average 



a 



(a^)hvp 



(694.9 ± 6.4) ■ 10^ 



-10 



(16) 
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There are some 0{a^) contributions already included in (0), namely, the 
intermediate 7r°7* and rj'y* states in a{e~^e~ hadrons). The rest of the 0{a^) 
corrections have been calculated in 1133, 13711 and are well under control 



a 



(q^) hvp 



;-io.o±o.6) ■ 10 



-10 



(17) 



It was pointed out in 
up to 80 % in the Novosibirsk data analysis of a{e~^e 
were estimated in and added back to their final number (|1T 



that final state radiative corrections are eliminated 

7* — > hadrons). They 
^ These 0{a^) 

corrections depend obviously on the experimental set-up and on the particular 
analysis of the data. In fact, they were already taken into account in the ALEPH 
tau dataQ used by and p8| and therefore should not be added back to the 
results (0) and (|I|). 

Summing (|1^) and (]T7|), one gets for the total hadronic vacuum polarization 
contribution 



hvp 



(684.9 ± 6.4) ■ 10" 



-10 



The expected e~^e~ data accuracy below 1 % supplemented with theoretical efforts 
like ||2^, ^] will reduce the uncertainty in from e~^e~ up to the order of 
6 ■ 10~^°. Joint works of tau and e'^e" groups ^ are also announced and 
will reduce this uncertainty further. Isospin breaking studies like will help 
to take these corrections under better control. Chiral symmetry can also help to 
reinforce the accuracy of the tttt dominant contribution [R^. 



4 Results and Summary 

Summing all the Standard Model contributions, (^), (^), (|]), and (|T^) to muon 
g — 2, one gets 

aj^ = (11 659 179.2 ±9.4) ■ 10"^°. (19) 

where uncertainties of the QED, EW, and hadronic light-by-light contributions 
have been added linearly and afterwards added quadratically to the hadronic 
vacuum polarization uncertainty. 
As a final result, we get 

a^- af^ = {23.1 ±16.9) -10-^^ (20) 

i.e. there is at present a bit more than one sigma of discrepancy. The signifi- 
cance of this discrepancy could be largely enhanced by the aimed experimental 
uncertainty BNL goal of 4 ■ 10~^°. The announced improvements on the hadronic 
contributions are also very interesting and can reduce the uncertainty of the muon 
g ~ 2 Standard Model prediction to the order of 6 ■ 10~^° to 7 • 10~^°. The near 
future of muon g — 2 reveals thus very exciting. 

thank Michel Davier for informing mc on this point. 
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